A possibility to measure elastic photon— photon scattering in vacuum 
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Photon-photon scattering in vacuum due to the interaction with virtual electron-positron pairs is 
a consequence of quantum electrodynamics. A way for detecting this phenomenon has been devised 
based on interacting modes generated in microwave waveguides or cavities [G. Brodin, M. Marklund 
and L. Stenflo, Phys. Rev. Lett. 87 171801 (2001)]. Here we materialize these ideas, suggest a 
concrete cavity geometry, make quantitative estimates and propose experimental details. It is found 
that detection of photon-photon scattering can be within the reach of present day technology. 
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I. INTRODUCTION 

Classically, photon-photon scattering does not take 
place in vacuum. However, according to quantum elec- 
trodynamics (QED), such a process may occur owing to 
the interaction with virtual electron-positron pairs. An 
effective field theory containing only the electromagnetic 
fields can be formulated in terms of the Heisenberg-Euler 
Lagrangian P, Q] , which is valid for field strengths below 
the QED critical field lO^^Vm"^ and for wavelengths 
shorter than the Compton wavelength 10~^^m 
Several suggestions to detect photon-photon scattering 
in laboratories have been made 0, 0) 13 j and the recent 
increase in available laser intensities has stimulated var- 
ious schemes ISl (l^, El 111 Ellll- In Ref. |3 we 
suggested an alternative method, using the significantly 
weaker (but still strong) fields that can be confined in 
a microwave cavity. The main advantage of adopting a 
cavity is that we can achieve a resonant interaction be- 
tween the eigenmodes in a large volume, leading to an 
excited signal at a new eigenfrequency. 

In the present paper we materialize the proposal made 
in Ref. [l^ in several ways. Various concrete geometries 
fulfilling all resonance and frequency matching conditions 
in a cavity are devised, and the coupling between the 
pump modes and the excited mode is evaluated for a rect- 
angular prism and a cylindrical geometry. The amplitude 
of the excited mode is then determined in terms of the 
quality factor of the cavity and the field strengths of the 
pump modes. Comparing with the performance reached 
in existing superconducting niobium cavities, it turns out 
that the cavity key parameters, namely the quality fac- 
tor and the allowed field strength (before field emission 
and/or superconductivity break down) can reach values 
which allow for several photons in the excited mode. To 
be able to detect the very weak excited signal in the pres- 
ence of the pump waves, we propose a "cavity filtering 
geometry" . Finite element calculations are made in or- 



der to demonstrate a geometry that fulfills the relevant 
resonance and frequency matching conditions. The esti- 
mated signal in the filtered region corresponds to roughly 
30 microwave photons/m"^. The results of Refs. and 
|l7j , where effects of single microwave photons confined in 
cavities are measured using sensitive techniques involving 
transitions of Rydberg atoms, suggest that detection of 
the estimated signal is possible. Hence we suggest that 
elastic photon-photon scattering can be observed with 
present day technology. 



II. BASIC EQUATIONS AND PRINCIPLES OF 
CALCULATION 



Photon-photon scattering due to the interaction with 
virtual electron positron pairs can be described by the 
Heisenberg-Euler Lagrangian [H 



i = eoF + keI [4F2 + 7G2] 



(1) 



where F — (E"^ 



c^B^)/2 and 
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G = cE • B. Here 
1.63 X 10"^" ms^/kg, a is the fine- 
structure constant, Ti Planck's constant, me the electron 
mass, and c the velocity of light in vacuum. The last 
terms in JQ) represent the effects of vacuum polarization 
and magnetization. The QED corrected Maxwell's vac- 
uum equations can then be written in their classical form 
using D = £oE + P and H = B/^o - M where P and M 
are of third order in the field amplitudes E and B, and 
jjLQ — 1/c^eo- Expressions for P and M can be found in, 
for example, Refs. 3, 4, 5, 6,|3,i,i,liol|l3,|r2, 13, 15]. 

There are several reasons for considering wave interac- 
tions in cavities: 

1. We can benefit from coherent resonant interactions. 
In contrast, the nonlinear coupling vanishes for par- 
allell plane waves in an unbounded medium. How- 
ever, we note that the presence of an inhomoge- 
neous background magnetic field can also be re- 
sponsible for a nonzero effect, see Ref. 01 • 

2. The growth of the new mode will not be saturated 
by convection out of the interaction region. 
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3. The techniques for detecting small signals in such 
cavities are very well developed, see e.g. Refs [a 

m 

Calculations of the coupling strength between various 
eigenmodes can be made including the nonlinear polar- 
ization and magnetization, see e.g. Ref. |l5j| . However, a 
more convenient approach, which gives the same result, 
starts directly with the Lagrangian density 

Our general procedure for finding the cavity eigen- 
mode coupling and the saturated amplitude of the ex- 
cited mode, to be applied in sections III and IV, can be 
summarized as follows: 

1. We determine the linear eigenmodes of the cavity 
using the standard Maxwell vacuum equations together 
with the boundary conditions for infinitely conducting 
walls, and express all field components in terms of the 
vector potential amplitude. 

2. Next we choose resonant eigenmodes fulfilling fre- 
quency matching conditions. For two interacting initial 
pump modes (indices 1 and 2) with distinct frequencies 
LUi and UJ2, the possible frequency matching choices, cor- 
responding to a cubic nonlinearity, are ±(2a;i ± 0^2) and 
±(wi ± 2^2). For definiteness we concentrate on the 
choice 



073 = 2ll!i — LU2 



(2) 



in all examples. Since the eigenfrequencies are deter- 
mined by the geometry, we note that the condition ^ 
gives a design requirement involving the dimensions of 
the cavity. In this context it can be noted that there are 
several reasons for using two pump waves rather than 
one, that in principle could excite a mode with three 
times the frequency of the pump mode. Firstly, for two 
pump modes there are generally better possibilities to 
vary the parameters in the set-up in order to optimize 
the output of the excited mode. Secondly, for a cylin- 
drical geometry and a single pump mode, the design re- 
quirements following from the frequency matching results 
in eigenmodes that have a too weak nonlinear coupling. 
Finally, there is a general tendency to get stronger non- 
linear coupling when the pump modes and the excited 
mode have rather close frequencies. 

3. We then perform a variation of the amplitude of the 
eigenmodes and let <5/ L (Prdt = in order to obtain the 
mode-coupling equations directly from the Lagrangian 
density 1^. The evolution equation for each eigenmode 
is obtained by expressing the Lagrangian in terms of the 
potential and varying the corresponding vector potential 
amplitude. The lowest order linear terms then vanish, 
since the dispersion relation of each mode is fulfilled. For 
the terms that are quadratic in the fields we must thus 
take into account that the amplitude has a weak time 
dependence when making the amplitude variation. How- 
ever, for the QED correction terms in the Lagrangian the 
time dependence of the amplitudes can be neglected. 

4. In the absence of dissipation, the equations now ob- 
tained imply steady growth of mode 3, until the energy 



of that mode is comparable to that of the pump modes. 
However, when some damping mechanism is present (e.g. 
due to a finite conductivity of the cavity walls) the ampli- 
tude saturates at a level where the mode-coupling growth 
balances the dissipation of the excited mode. This ef- 
fect can be easily included by adding a phenomenologi- 
cal damping term in the evolution equation, whose value 
is estimated by comparing with quality factors currently 
reached in superconducting niobium cavities. 



III. RECTANGULAR PRISM GEOMETRY 

We start by considering a rectangular prism cavity, 
with one of its corners in the origin, and the opposite 
corner with coordinates {xo^yo, zq). In practice we are 
interested in a shape where zq » xq, yo but this assump- 
tion will not be used in the calculations. We let the large 
amplitude pump modes have vector potentials of the form 



Ai = ^1 sm I — sm exp(— icjiijy + c.c (3) 

XoJ \ Zq 



and 



A2 = A2 sm — sm exp(— ia-'2i)x -f c.c. (4) 

\VoJ \ zo J 

where c.c. denotes complex conjugate, ni_2 = 1,2,3..., 
and where we have chosen the radiation gauge such that 
the scalar potential is zero. It is easily checked that the 
corresponding fields (omitting the c.c.) 

Biz = — Ai COS — sm exp(— ia;i<), (5a) 

xo \xoJ \ zo J 

„ niTT / TTX\ ( n\'Kz\ 

Bix = Ai sm — cos exp(— itJitX,5b) 

Zq \xoJ \ Zq J 

f ■nx\ ( n\-Kz\ 

Ely — luJiAi sm — sm exp(— iwit), (5c) 

\xoJ \ ZQ J 

together with ujf = n^-K^ (? j zf^ + tt^c^/x^, and 

' 1^y\ . ( 71.2 TTZ 



B2z = A2 cos — sin 

yo Kvq ) \ Zq 



exp(— ia;2i))(6a) 



naTT . . / Try \ / n2'Kz , , ■ iau\ 
B2y = A2 sm — cos exp(-iW2i),(6b) 

Zo \yoJ \ Zq ' 





' n27rz\ 








K ZQ J 



exp(-iw2<), (6c) 



together with aj| = ti^tt^c^ /zq + tt^c^/j/q are proper 
eigenmodes fulfilling Maxwells equations and the stan- 
dard boundary conditions. Similarly we assume that the 
mode to be excited can be described by a vector potential 



A A ■ ( ■ ( "-37rz\ 
A3 = A3 sm — sm exp 

\VqJ \ zo J 



)(— ia;3t)x + c.c. (7) 
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where a;| = n^-K^c^ / Zq + tt^c^/j/q, in which case we get 
fields of the same form as in Eqs. (|6a|) - (|^ . 

Next we turn to point three in the scheme of the previ- 
ous section. As noted above, when performing the varia- 
tions 5A'^, the lowest order terms proportional to 5A^A^ 
vanish due to the dispersion relation, and we need to 
include terms due to the time dependence of the am- 
plitude of the type A'id{5A%^)/dt. For the fourth order 
QED corrections proportional to 5A^, only terms pro- 



portional to A^AjMg survive the time integration, due 
to the frequency matching ||2Jl. After some algebra the 
corresponding evolution equation for mode 3 reduces to 



dAs _ _i£oKw| ^ .2 4* 
dt ~ 8 ^-^^1^2 



where the dimensionless coupling coefficient K-^cc is 



(8) 




8^ 

2 2 
^0% 



-)uj2ni,{- 



Tni 



V 



CJ2'^3 



r^2?^37^'c'' / ln\ 



(9) 



The three different sign alternatives in © correspond to 
the mode number matching options 



2ni — 712 
2ni -t- n-2, 
2ni — 712 



ns = 
ns = 
ns = 



(10a) 
(10b) 
(10c) 



respectively, which must be fulfilled in order for the cou- 
pling to be nonzero. It is now possible to evaluate the 
coupling coefficient for specific mode numbers and ge- 
ometries consistent with the frequency matching condi- 
tions. As described in the previous section, we could 
then determine the saturated amplitude from a balance 
between mode-coupling growth and dissipation, where 
the latter effect can be introduced by simply adding a 
phenomenological damping term. However, it turns out 
that a cylindrical geometry gives a slightly better perfor- 
mance, and thus we will instead work out that case in 
more detail. 



IV. CYLINDRICAL GEOMETRY 

We consider the case where all eigenmodes are TE- 
modes with no angular dependence, with fields that can 
be derived from the vector potential 



where a is the cylinder radius, zq the length of the cavity, 
Ji the first order Bessel function and (3 one of its zeros. 
The cylinder occupies the region < z < zq centered 
around the z-axis. We have here introduced cylindrical 
coordinates p and z as well as the unit vector ip in the 
azimuthal direction. The corresponding fields are 



E = iujAJi{pl3/a) sin ( 1 exp(-iwt)^ -I- c.c(12) 

PA , , „ I s . ( n-Kz \ , . , ^ 

B = Jo(pp/a)sm exp(— i(jji)z — 

a \ J 

J17T ( T17T Z \ 

— AJi{pP/a) cos exp(— iLjt)p + c.c. (13) 

zo \ J 



where the eigenfrequency is uj^ — c^[(/?/a)2 + (mr/zo)'^] 
for all modes uj = 0^1,2,3. We note from the frequency 
matching condition lut, — 2lui — 1^2 that all of the eigen- 
modes cannot have the same order of their respective /3, 
for resonant interaction to occur, and thus we introduce 
(3 — /9i.2,3 for the different modes. Proceeding in the 
same manner as described in the two previous sections, 
we obtain after lengthy but straightforward algebra 



dAs 
dt 



ieoKwf 



Kcy\A^A2 



(14) 



A = AJi{pf3/a) sin 



TITTZ 



exp(— itjt)^ -f c.c. (11) where the cylindrical coupling coefficient Kcy\ is 



8c2 



cyl 



27r 71^712713 



a + {I3in2m{la + Ic) + 2/327ll7l3/d - 2/3371i7l2/e) 



a^z'^ 



{'^P'iU}iUJ2le + 2/32WlW3/<j - /3iW2CJ3(/c + 34)) 



(?l3WiCi;2 - «2^^lW3 - '^l'^2^^3) 



(15) 



4 



and the integrals are defined as 



la - / Jt (M Ji {P2U) Ji iPsu) u du, (16a) 
Jo 



h = / [3 J2 (/3im) + J? 
Jo 

X Jo i02u) Jo iPau) u du, (16b) 



-fc = / Jq ifiiu) Ji{P2u) Ji{l33u)udu, (16c) 

JO 



Id= / Jo(/3i?i) Ji(/3iu) Jo(/32u) Ji(/33M)udw, (16d) 
Jo 

and 



4 = / Jo iPiu) Ji {f3iu) Ji{P2u) Jq {P3u)udu. (16e) 
Jo 

When calculating H15() we have assumed the mode num- 
ber matching ^3 = 2ni + n2. Equation (|14|l implies a 
linear growth of mode 3, until the backreaction of the 
pump modes becomes significant. In reality dissipative 
mechanisms (e.g. a finite conductivity of the cavity walls) 
will put a limit on the excited amplitude. This can 
be described in a phenomenological way by substituting 
d/dt — > d/dt — {uj3/2'kQ), where Q is the cavity quality 
factor. The steady state amplitude is thus 



char 



(17) 



where we have introduced the characteristic electric field 
£:^har = (eok)"^/^ « 2.6 X lO^Oy/m [l^. Looking at 
the number of excited photons in the cavity mode N « 
£0 J EaEsd'^f /fi^3 we finally obtain 



N, 



(18) 



Before evaluating IjlSfl we need to specify the mode 
numbers and the geometry. As an example we let 
(ni, 712,^3) — (3,15,21) (fulfilling — 2ni + 71,2) and 
take ^2=^3 — 3.83, corresponding to the first zero of 
Ji, and /3i = 7.01 corresponding to the second zero. This 
gives us zo/a — 9.53 through the frequency matching 
condition ^ and determines the frequency relations to 
'-^3/^2 = 1-26 and UJ3/UJ1 = 1.12. Substituting these val- 
ues and numerically evaluating the integrals l|16a|) - (|16e|) 
then gives -f^cyi = 3.39. The key parameters are the qual- 
ity factor and the pump field strength. An advantage 
with our choice of eigenmodes is that the pump electric 
field is zero at the cavity surface, which means that we do 
not have to worry about field emission Instead the 



pump amplitude is limited by the surface magnetic field, 
which needs to be below the critical value for which the 
walls cease to be superconducting. From the experimen- 
tal results presented in Ref. [23| , we find that the critical 
magnetic field for a high pure niobium material can reach 
B « 0.28 Tesla, at a temperature of the order of or below 
1 K. If we specify zq = 2.50 m, we get uii = 8.10 x 10^ 
rad/s and ^2 = 7.15 x 10^ rad/s, implying that field levels 
close to the magnetic surface field condition correspond 
to Ai = 0.017 Vs/m, and A2 « 0.024 Vs/m. The con- 
ductivity of niobium allow for quality factors Q > 10^^ 
for temperatures in the range of interest below 1 K. But 
we note that these high levels have been hard to reach 
in practice, although Q ~ 4 x 10^° in Ref. 0. We also 
note that there has been a tendency to get lower quality 
factors when applying stronger fields. However, in Ref. 
|2ll | , it has been shown that it is possible to reach surface 
fields of the order of the critical level without significant 
decrease of the Q-factor. Thus adopting Q = 4 x 10^° 
and the rest of the parameter values as specified above, 
we obtain 



QED 



18 



(19) 



Keeping the cavity at a temperature T w 0.5 K, thus 
means that the number of generated photons in our ex- 
ample is well above the thermal fluctuation level TVth = 
l/[exp(?iw3/fcT) — 1] R:! 7, where k is the Boltzmann con- 
stant. To get an even lower thermal fluctuation level, it 
would be of interest to generate photons with higher fre- 
quencies. However, it is not wise to just scale down the 
dimensions used above to get Hivs/kT > 1 and A'^th ^ 1, 
since A^qed decreases too quickly with the cavity volume. 
Instead we could consider higher mode numbers '1.1,2,3 
and higher orders of /3 to get a larger excited eigenfre- 
quency fulfilling huJs/kT > 1 which gives a higher ratio 

NQED/Nth. 

On the other hand, detecting a very weak signal in the 
presence of strong pump fields might be difficult even 
for signals well above the thermal fluctuation level, and 
thus we will below investigate the possibility of a cavity 
geometry that directly filters away the pump signals. 



V. CAVITY FILTERING GEOMETRY 

We search for a cavity geometry with the following 
properties: 

1. The cavity should be rotationally symmetric 

2. There should be three eigenmodes with frequencies 
fulfilling 

3. There should be a cylindrical interaction region 
where the mode structure resembles that consid- 
ered in the previous section. 

4. There should be a filtering region with a cross sec- 
tion small enough to give an exponential decay of 
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the pump modes, but large enough for the excited 
eigemnode to propagate. 

5. There should be an entrance region for the pump 
modes. 

We thus make a geometrical design of the type out- 
lined above, see Figs. la,b,c. The eigenmodes are then 
calculated using the method of finite elements. In gen- 
eral there will be eigenmodes with the desired proper- 
ties, except that they will not fulfil Eq. (0) for the spec- 
ified geometry. However, by repeated calculations vary- 
ing the length I of the cavity, the mismatch frequency 
Suj{l) = uj3{l) — 2uji{l) +ll>2{1) gradually approaches zero. 
The dimensions of the cavity and the corresponding mode 
structure is shown in Figs 1 a,b,c. As can be seen, the 
pump modes have decreased their amplitudes by a factor 
of the order of 10®, from the interaction region to the 
end of the filtering region (to the left), whereas mode 3 
has essentially the same amplitude in both regions. By 
increasing the filtering distance, the pump signals could 
of course be reduced further. Naturally a filtering ge- 
ometry will make the mode coupling somewhat weaker, 
as compared to our pure cylinder example. By adding a 
filtering region of roughly the same size as the coupling 
region, the coupfing factor reduces to Kfn ~ Kcyi/V^, 
in which case we keep the same number of excited pho- 
tons if we choose the size of the total cavity region to be 
roughly twice the size of the cylinder cavity presented in 
the previous section. 



VI. NONLINEARITIES IN THE WALLS OF 
THE CAVITY 

To our knowledge, a well-established and simple non- 
linear model for the superconducting RF-state does not 
exist. As a starting point for a discussion, we may adopt 
a model with a nonlinear magnetic third order suscepti- 
bility AP = x^^'"'^'' BjBkBi, using the Einstein summation 
convention, and defining the susceptibility in terms of B 
rather than H. We then consider the same geometry of 
the cavity and the eigenmodes as in section IV. The mag- 
netic z-components of the pump fields penetrate roughly 
a skin depth inside the walls. For a nonzero value of 
^3333 = ^^j^ ^Yie part of the nonlinear magnetization that 
can act as a source for mode 3 is then M3 = Xn\Bi^B2zZ. 
Acoordingly we get currents in the azimuthal direction 
Jsip'fi = V X M3, which in turn may induce radial vari- 
ations in the magnetic field B^^ip) = (mo/p) / J^ippdp. 
At the same time, the jump in the nonlinear susceptibil- 
ity across the vacuum/superconductor boundary causes 
a jump in the magnetization, and thereby a surface cur- 
rrent jas — — M3 x p. A combination of the surface and 
bulk currents then yields a nonlinearly induced magnetic 
field inside the vacuum region 



^3zL<a = (/"O/P) / JsvPdp- fJ-OXnlBl^B. 



2z\ 



(20) 



where A should be much larger than the skin depth such 
that JstpiX) is negligible. Carrying out the integration 
in Eq. (|20|) . we see that the two terms cancel such that 
B^zlp^a ~ i-^- nonlinear currents do not induce a 
magnetic field inside the cavity. 

We stress that this is not a proof that nonlinearities 
in the walls generally are unable to affect the physics 
inside the cavity. However, it suggests that a nonlinear 
current in the walls with the proper eigenfrequency does 
not by necessity excite the corresponding eigenmode of 
the cavity. 



VII. DISCUSSION AND SUMMARY 

In the present paper we have materialized our proposal 
for the detection of elastic photon- photon scattering . 
In particular, we have calculated the output level of scat- 
tered photons in terms of the allowed pump field strength 
and the cavity quality factors for a reactangular prism as 
well as for a cylindrical geometry. Furthermore, we have 
made finite element calculations to show that the reso- 
nance and frequency matching conditions can be fulfilled 
in a filtering geometry, where only the scattered mode has 
a high enough frequency to reach the cavity region with a 
lower cross-section. Using performance data from current 
state-of-the-art superconducting niobium cavities, where 
a high quality factor Q w 4 x 10^" is combined with sur- 
face field strengths close to the critical value B « 0.28 
Tesla, we deduce that the number of scattered photons 
can reach A^ « 18 in a cylindrical cavity with length 2.5 
m and 25 cm radius. In a filtering geometry we estimate 
that the number of photons/ volume will be reduced by a 
factor of order 2, as compared to the pure cylinder case. 
Recent results [l^ IT^ . where measurements involving 
transitions in Rydberg atoms interacting with single mi- 
crowave photons have been made, strongly suggest that 
the estimated field levels are within the range detectable 
with present day technology. 

We note that, in principle, nonlinearities in the walls 
of the cavity may lead to excitation of the same mode 
as caused by the QED nonlinearities. Our simple model 
calculation in section VI suggests that the mode-coupling 
due to such an effect vanishes. However, a more rigorous 
treatment is necessary to draw definite conclusions. QED 
theory predicts a definite output level of mode 3. Ex- 
periments that result in a much higher level would thus 
indicate that nonlinearities in the walls play the main 
role. 

Finite element analysis, may provide a starting point 
for a suitable design of the cavity. However, the degree 
of fine tuning of the resonance frequencies of the modes 
is very high, since for optimal performance the mismatch 
of the eigenfrequencies should not exceed Slo ~ lo/Q. 
Hence the final adjustments of the cavity geometry must 
be made experimentally. 
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FIG. 1: The geometry of the cavity and the mode structure for the filtering geometry. Note that only half of the cavity is 
shown, since the other half is redundant due to rotational symmetry of the cavity as well as the fields. The small region to 
the right is the entrance region, the large middle region is where the interaction takes place, and the region to the left is the 
filtering region. All modes have the electric field in the angular direction, i.e. E = E{p, z)ip. The variations of In \E{p, z)\ are 
shown in color code. 

a) The mode structure of pump mode 1. The exponential decay in the region of small cross-section diminishes the amplitude 
by a factor 10~^ in the end of the filtered region. 

b) The mode structure of pump mode 2. The exponential decay in the region of small cross-section diminishes the amplitude 
by a factor 10~* in the end of the filtered region. 

c) The mode structure of the excited mode. The amplitude is roughly the same in the region of interaction and in the filtered 
region. 



